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Abstract Gluon propagator in the infrared limit, computed in the Landau gauge, 
can be cast into an universal form with a running mass for the gluon. In this 
way we are able to show that all the proposals appeared so far in literature are 
equivalent and describe the same physics. In this way, it appears essential to reach 
a general agreement about the interpretation of such results. We discuss the points 
that make difficult to agree on a similar view about this. 


1 Introduction 

Since the initial proposal by Gribov and Zwanzinger mm, the study of correlation 
functions, mostly the propagators, of the Yang-Mills theory in the Landau gauge 
has seen a lot of activity both theoretical and through lattice computations. Gribov 
and Zwanzinger yielded arguments favoring a gluon propagator running to zero 
and a ghost propagator running to infinity faster than a free one when momenta 
go to zero, the infrared limit. This proposal was crucial to understand confinement 
^e-mail: marcofrasca@mclink.it 






2 


as, being this the state of art of our understanding, all the criteria for a conhned 
Yang-Mills theory were promptly satisfied. Indeed, this proposal entails a non 
trivial infrared hxed point for a pure Yang-Mills theory in absence of quarks. Few 
years later, this area of research has known a signihcant growing interest thanks to 
the work by Alkofer, von Smekal and Hauck [3113] . that used a truncated scheme to 
solve Dyson-Schwinger equations and by Cornwall [5] , that introduced the so-called 
pinch technique in perturbation computations [6] . But while the work by Alkofer, 
von Smekal and Hauck supported the conclusions by Gribov and Zwanzinger, 
Cornwall’s approach obtained a significant different result: Gluon propagator did 
not seem to go zero rather it reached a finite non-null value in the infrared limit. At 
that time, lattice computations were not able enough to decide what was the right 
behavior. On the other side, some other theoretical works supported Cornwall’s 
conclusions [TlIHllQlfTn] following different approaches. 

Finally, on 2007, due to the improvement of computational resources, it was 
possible to reach volumes large enough to verify the behavior of propagators in 
the Landau gauge at lower momenta doing lattice computations. The results con¬ 
firmed that the gluon propagator reaches a finite non-null value at lower momenta 
and the ghost propagator behaves like that of a free massless particle miMM- 
Thereafter, the original proposal by Gribov and Zwanzinger was properly corrected 
to £t lattice data m as also confirmed from lattice computations m- Today, this 
result has been widely accepted for three and four dimensions while for the two- 
dimensional case the original scenario by Gribov and Zwanzinger, dubbed scaling 
solution, applies [16]. This means that the scaling solution is not meaningful for 
our analysis but could have a meaning in three and four dimensions in QGD when 
quarks are present as recently shown in m- 
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Our aim in this paper is to show how all the theoretical approaches presented 
in literature so far to analyze propagators in Landau gauge for Yang-Mills the¬ 
ory without quarks yield identical results in the deep infrared. This means that, 
whatever approach one likes, it is possible to do any kind of computations with 
the theories devised so far. In order to achieve this conclusion, we just use the 
original idea of Cornwall [5] that, in the infrared limit, the gluon propagator is 
characterized by a running mass that reaches a constant value at very small values 
of momenta. This idea applies straightforwardly to whatever approach was con¬ 
ceived so far to discuss Yang-Mills theory in the infrared limit. The only critical 
point is to understand the meaning of such a result in terms of running coupling. 
There are diverging views about this but the facts that perturbation theory ap¬ 
plies successfully in the infrared limit, as shown by Tissier and Wschebor nmg, 
and renormalization group arguments confirm the triviality of the theory in the 
infrared limit, as shown by Weber [20], seem to seriously hint that the running 
coupling should go to zero lowering momenta. This was seen on the lattice m 
with the definition of the running coupling given therein. 


The paper is so structured. In Sec. [2] we introduce the gluon propagator in the 
Landau gauge and show how the different approaches yield it in closed analytical 
form. In Sec. |3| we give a dehnition of running mass and evaluate it for all the 
approaches we discuss. The equivalence is seen at numerical level between them 
being perfectly coincident below 1 GeV. In Sec. Hjwe discuss the running coupling 
and yields a triviality argument depending on the preceding analysis of the running 
mass. Finally, in Sec. [5] conclusions are presented. 
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2 Gluon propagator 


In the Landau gauge, the gluon propagator takes the generic form 



( 1 ) 


being a, b color indexes and the metric tensor. In the proposals seen in current 
literature, the gluon propagator in the Landau gauge can be generally represented 
through the equation 


Zl ^{p) = Z + 


( 2 ) 


where we assume an Euclidean signature, M(p) is a running mass and Z a constant 
factor. The running mass has the property to reach a finite value when p —>• 0, in 
agreement with lattice data, and becomes negligibly small in the opposite limit 
|22| . This formula yields also a definition of a running mass by inversion as 


(p) = -p^ + ZA ^ (p) 


(3) 


that can be immediately used in numerical calculations. The evaluation of the 
running mass is meaningful just below about 1 GeV where the coupling is large 
enough that perturbation theory cannot be trusted anymore. In all the studies 
given in literature the constant Z can be taken to be unity. The running mass will 
depend on the gauge choice as it depends on the gluon propagator that changes 
with the given gauge. 

The current proposals presented in literature, by different approaches, can be 
resumed as follows. The propagator proposed by Cornwall [5] and then refined in 
a series of works by Papavassiliou, Aguilar and Binosi can be cast in the form |231 
[ 21 ] 
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where m, g and p are taken as free fitting parameters and p is the renormalization 
point taken to be 3 GeV. From this form of the propagator the derivation of the 
running mass is straightforward and is given by 


^CPAb{p) — rri^ +P^ 


9671^ 


■In 


2 I 2 

p + pm 


(5) 


We assume this formula to hold up to p « 1 GeV. This expression holds in the 
Landau gauge and will change with a different choice of the gauge [23] . 

In the infrared limit, Tissier and Wschebor used perturbation theory after 
inserting a mass term for the gluon in the Lagrangian. Their propagator takes the 
form [TSlIT^ 


. 2 2 p^N 

^Tw(p) = P + QOA- 
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( 6 ) 


and the running mass can be immediately obtained as 


Mrwip) = + 


P^Ng^ 
3847r2 
2 N 3 
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(7) 


This equation was obtained from perturbation theory and, as such, could be im¬ 
proved by higher order corrections. The parameter m is the gluon mass introduced 
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initially in the Lagrangian and should be fixed to fit while ^ is the renormalization 
point. 

Postulating the existence of a dimension-two condensate, Dudal, Gracey, Sorella, 
Vandersickel and Verschelde proposed a gluon propagator in the form m 

A / \ ^0 /o\ 

DGSVV P _|_ 2Ng'^'y‘^ + m'^M^ 

being mo, Mq and 7 free fitting parameters. A study has been accomplished com¬ 
paring this form with lattice data and the best fit was obtained for 


^DGSVvip) 


0‘+ , 


( 9 ) 


that can be easily rewritten in to the form given in eq.®. But this is also a 
very good approximation to the propagator we proposed some time ago m and 
recently refined [26] in the form 


Af{p) 


00 


E 

n=0 


Bn 


( 10 ) 


where we put 


Bn = i2n+iy 




4A'3(-1) 1 -pe-( 2 "+i)’r- 


( 11 ) 


and rrin = (2n + 1) 2 x^- 1 ) H®''® ™-o is a free fitting parameter and K{—1) is 
the complete elliptic integral of the first kind. It is easy to see that each term 
in (nni) is exponentially damped and so, the propagator in eq.(|9|) is an excellent 
approximation to it. In the following we will assume that the parameters of eq. © 
are always chosen to best fit eq. (nni). The reason for doing so is that, based on the 
preceding observation, these two forms of propagator should give the same results 
when used to fit lattice data and the running mass is almost the same in the 
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momenta region we aim to stndy. In any case, the rnnning mass can be evalnated 
using eq. also noticing that 

OO 

Y^Bn = l. (12) 

n=0 

This propagator is represented in Fig. [T]and is similar to that observed on lattice 
computations that also all other proposals excellently fit. 



Fig. 1 Numerical evaluation of the gluon propagator in the Landau gauge as given in eg. lllOII . 
for rriQ = 0.7 GeV, in close agreement with lattice computations. 


3 Running mass 

When evaluated, running masses for all the devised approaches yield the cnmula- 
tive behavior displayed in Fig. [2l Here TW means Tissier and Wschebor, CPAB 





means Cornwall, Papavassiliou, Agnilar and Binosi and DC means Dudal and Cuc- 
chieri referring to the numerical fits in for the propagator of Dudal, Gracey, 
Sorella, Vandersickel and Verschelde. 



Fig. 2 Numerical evaluation of the running masses for all the approaches discussed in Sec. [5] 
(see text for further details). 


In Fig. [3] we show just the fit with the CPAB running mass of the running 
mass given by eq. (HOI and it is clearly evident the excellent agreement in the 


region below 1 GeV. 
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Fig. 3 Numerical evaluation of the running masses just for CPAB and from propagator lllOII . 


These fits are provided with the parameter m in eq. dZD and eq. (ED to be the 
same as mo in eq. (nnD. In any case we assume for the renormalization point = 
3 GeV. This applies both to TW and CPAB. For CPAB we also take p = 27 and 
Qq = 9.77. We assume SU(3) everywhere so that N = 3. All the running masses 
agree excellently till momenta around 1 GeV as expected except for TW. Latter 
discrepancy is due to the fact that TW running mass is obtained by perturbation 
theory at one loop and so is amenable to improvements going to higher orders. 
We expect in the end an excellent agreement for all the approaches discussed here. 
This result has some implications. Firstly, whatever the preferred approach the 


results are always in agreement with lattice data but in some cases there are more 








10 


parameters to adjust the fit. This is not relevant in view of the fact that the 
qualitative picture they give is a correct one. Secondly, the running mass given in 
CPAB, eq.Q, is a really good fit for all the cases we discussed so far and represents 
an excellent resummation for eq. (nni). Thirdly, it appears a fundamental fact that 
a better interpretation of this equation for the running mass will give a better 
understanding of Yang-Mills theory in the infrared. The key point here is to reach 
an agreement about the theory being or not trivial in the infrared limit. Common 
wisdom assumes that the theory should reach a finite non-null infrared fixed point 
and the coupling should have a finite non-null coupling at zero momenta. There 
is no proof whatsoever of such a statement for a pure Yang-Mills theory while 
this is certainly true for QCD when also quarks are included. The main difficulty 
to proof such a statement relies on the difficulties on share a common view on 
how an infrared running coupling should be properly defined. Fourthly, all these 
approaches draw the important conclusion that Yang-Mills theory has indeed a 
mass gap but none of them provides a really complete mathematical proof. 


4 Running coupling 

There are some important results that add to the conclusion that Yang-Mills 
theory becomes free in the low-energy limit without quarks. We try to summarize 
them here. As discussed in the preceding section, Yang-Mills theory with a mass 
term added reproduces quite well lattice data using standard perturbation theory 
nails]. Recently, it was shown by Weber and Dall’Olio m that this perturbation 
technique can be used to solve Callan-Szymanzik equations for the theory in the 


infrared. The running coupling is shown to go to zero at lower momenta as already 
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shown by Weber using renormalization group arguments [20]. This would be in 
perfect agreement with lattice studies (see m and Refs, therein) and the rather 
commonly accepted definition of the running coupling as 

Mp) = ^Z{p)j\p) (13) 

assuming for the gluon propagator A{p) = Z{p)/p^ and for the ghost propagator 
G{p) = J{p)/p^. It easy to see that, with this definition, all the cases we discussed 
in Sec. |2]will have the running coupling going to zero as momenta go to zero. This 
is due to the mass gap that all these approaches share and with the fact that, for 
the ghost propagator, J{p) « 1 when p ^ 0. It is also fundamental to note that 
the running mass has the property that 

M‘^{p)^^rTn? (14) 


and so we are left with a Yukawa propagator proper to a free theory in agreement 
with triviality requirements. So, the existence of a mass gap implies that the theory 
would be trivial in the infrared limit. 

Finally, let us consider a propagator in the form 




mil + ie 


(15) 


with Bn > 0 for any n and mn the mass spectrum. In this case, Kallen-Lehman 
representation can be applied with 


= '^ZnS{fi^-ml). (16) 

n=0 

Given the case that = 1, the theory has no bound states and so it is 

trivial. This is exactly the case provided in eq. (nni) that fits excellently well with 
the other proposals in the low energy limit in agreement with Ref. |28| for scalar 




12 


field theory. Indeed, we can write the approximate equation at small momenta, in 
the Euclidean case and using eqs. (Uni), (HH), the corresponding spectrum for rrin 
and eq. ©, 


2 I 

-P + 


E 

n=0 


Bn 


+ mi 


2^ 2UNg\ 
rn + p „„ In 


9671 ^ 


2 I 2 

p + pm 


(17) 


provided 


— I 

Vn=0 ” 


(18) 


Then, the ghost sector decouples from the the Yang-Mills field and a free massless 
propagator holds in the deep infrared limit. We can evaluate the running coupling 
through eg. (11311 by noticing that J{p) « 1 and 


Z(P) 


p-+ In 


(19) 


and then 



that goes to zero like p^ as seen on lattice m using the same definition of the 
running coupling. This fully justifies the successful perturbative approach of Tissier 
and Wschebor in the infrared limit mm- 

Quantum chromodynamics without quarks does not exist in nature. This means 
that a Yang-Mills theory could be trivial but adding quarks can change the char¬ 
acteristic of the infrared fixed point. This same thing is seen in the ultraviolet 
limit where the beta function depends critically on the number of quarks. For the 
latter case one has 


/3(as) 



( 21 ) 


where N is the number of colors and Nq is the number of quarks. We see that 
asymptotic freedom is preserved in QCD only if Nq < llN/2 otherwise quarks 
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could change the characteristics of the theory at higher energies. This is exactly 
what happens in the infrared limit as has been recently shown in Ref. EZl by 
solving Dyson-Schwinger equations. 


5 Conclusions 

The analysis performed for different gluon propagators proposed in literature so 
far has shown how these apparently different approaches should be considered 
equivalent below 1 GeV instead. The idea is that anyone of these could be refor¬ 
mulated through a proper definition of a running mass. This idea was firstly put 
forward by Cornwall on 1982 together with the view that the gluon propagator 
should display a mass gap when momenta are small enough. This scenario appears 
vindicated by lattice computations. Anyhow, it is interesting to note that the idea 
that a pure Yang-Mills theory should hit a non-trivial fixed point in the infrared 
limit, without any kind of proof whatsoever, was so rooted that, even if Cornwall 
obtained a Yukawa propagator in the deep infrared, no argument of triviality was 
put forward at all. Since then, arguments are mounting about this conclusion and, 
for some authors, this is no more a taboo. This will be material for sociology of 
science in the future. 

The most important conclusion to be drawn is that we have now a great arsenal 
of ammo to treat Yang-Mills theory in the infrared limit that is worth to be 
exploited and applied to more involved problems of quantum chromodynamics. 
My hope for the future is to see the low-energy spectrum of the theory completely 


under control. 
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